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Abstract 



Since the x dependence of the axial- anomaly effect in inclusive polarized deep inelastic 
CN . scattering is fixed, the transformation from the MS scheme to different factorization schemes 

cn ' are not arbitrary. If the quark spin distribution is demanded to be anomaly-free so that it 

fr^ I does not evolve with Q^ and hard gluons contribute to the first moment of gi{x), then all 

C^^ ■ 

t;;-!- . the moments of coefficient and splitting functions are fixed by perturbative QCD for a given 



75 prescription, contrary to the commonly used Adler-Bardeen (AB) or AB-like scheme. It 
is urged that, in order to correctly demonstrate the effect of factorization scheme depen- 
^ ■ dence, the QCD analysis of polarized structure functions in next-to-leading order should be 

I ' performed, besides the MS scheme, in the chiral- invariant factorization scheme in which the 

O-t. 

D , axial anomaly resides in the gluon coefficient function, instead of the less consistent and 



ambiguous AB scheme. 



Because of the availability of the two-loop polarized splitting functions APj • (x) recently 
|1]|, it became possible to embark on a full next-to- leading order (NLO) analysis of the 
experimental data of polarized structure functions by taking into account the measured 
X dependence of Q^ at each x bin. The NLO analyses have been performed in the MS 
scheme and the Adler-Bardeen (AB) scheme [2-10]. Of course, physical quantities such as the 
polarized structure function gi{x) are independent of choice of the factorization convention. 
Physically, the spin-dependent valence quark and gluon distributions should be the same in 
both factorization schemes. The recent analysis by the E154 Collaboration |§] has determined 
the first moments of parton spin densities in both schemes: 
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\'^"'v)mS ~ '-'•"^-0.02-0.04-0.01 ' y'^U'vJAB — "-'•'^-0.02-0.03-0.01 ' 
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where errors are statistic, systematical and theoretical, and AE = Au + Ad + As with 

Aq = Aq^ + Aqs = / dxAq{x) = / dx\q^{x) + q'{x) — q^{x) — q\x)\. (2) 

JO JQ 

We see that although the results of the fits in both MS and AB schemes are consistent 
within errors, the central values for Ag^, especially for AG, are not the same and the fits 
are significantly less stable in the AB scheme. This sounds somewhat annoying since if the 
polarized structure function is truly factorization scheme independent, then it is expected 
that the central values and errors of Aq^ and AG in the MS and AB prescriptions should be 
quite similar and that AS obeys the relation 

ASj^= ASab-?^AGab. (3) 

Also, because the NLO spin-dependent splitting functions AP/^ '{x) are originally calculated 
in the MS scheme, one may wonder if the polarized splitting kernels in NLO proposed in the 
AB scheme will render the evolution of gi{x, Q"^) scheme independent. 

In this short Letter, we wish to emphasize that since the x dependence of the axial- 
anomaly effect in the quark spin distribution or in the gluon coefficient function, depending 
on the chosen factorization scheme, is fixed, the transformation of coefficient and splitting 
functions from the MS scheme to the improved parton model scheme in which the spin- 
dependent quark distribution does not evolve with Q^ and hard gluons make contributions to 



the first moment oi gi{x), is also determined. As a consequence, we urge that a NLO analysis 
of (7i(x, Q^) data should be performed in the so-called chiral-invariant factorization scheme 
to be introduced below in order to see if gi{x, Q"^), Aq^^lx, Q^), AG{x, Q"^) are really scheme 
independent. Although none of the material presented in this Letter is new, a clarification on 
this issue is fundamentally important for the QCD analysis of polarized structure functions 
in NLO. 

Including QCD corrections to NLO, the generic expression for the polarized proton struc- 
ture function has the form 

9l{x,Q^) = ^Y.4[^C,ix,as)0Aq{x,Q')+ACGix,as)0AG{x,Q^) 

= ^Ee?[V°H^,Q')+A?^'n^,Q')+A?f(^,Q') (4) 

+ACJj^^ (x, as) ® Ag(°) (x, Q^) + ACg^ (x, a,) ® AG(x, Q^) + 

where uses of AC^°^(x) = 6(1 — x) and AC^^ (x) = have been made, ® denotes the 
convolution 

/(X)®^(X)= f'^f(^]g{y), (5) 

J^ y \yj 

and ACq, ACq are short-distance quark and gluon coefficient functions, respectively. More 
specifically, AC^^ arises from the hard part of the polarized photon-gluon cross section, 
while AC^^^ from the short-distance part of the photon-quark cross section. Contrary to the 
coefficient functions, Ag^(x) and Ag'^^^(x) come from the soft part of polarized photon-gluon 
and photon-quark scatterings. Explicitly, they are given by 

Ag(i)(x, Q') = A0(y^(x) ® Ag(°)(x, Q^), Agf (x, Q') = A0«e(x) ® AG(x, Q^), (6) 

where A0j/j(x) is the polarized distribution of parton j in parton i. Diagrammatically, A(p , 

q/G 



and A(j) L are depicted in Fig. 1 (see e.g. [|TT 



Since A^^^^ is ultravioletly divergent, it is clear that, just like the case of unpolarized 
deep inelastic scattering (DIS), the coefficient functions ACg and ACq depend on how the 
parton spin distributions Acjr-l are defined, or how the ultraviolet regulator is specified on 
A(l)^^\ That is, the ambiguities in defining Acp J and Acp L are refiected on the ambiguities 
in extracting AC^^ and AC^^ . Consequently, the decomposition of the photon-gluon and 
photon-quark cross sections into the hard and soft parts depends on the choice of the factor- 
ization scheme and the factorization scale /i [for simplicity, we have set /i^ = Q^ in Eq. (4)]. 
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Figure 1: Diagrams for the quark spin distributions inside the parton: Acj) , and Acj) ,q. 

Of course, the physical quantity g{{x) is independent of the factorization prescription (for a 
review on the issue of factorization, see |T2i). 



However, the situation for the polarized DIS case is more complicated: In addition to 
all the ambiguities that spin-averaged parton distributions have, the parton spin densities 
are subject to two extra ambiguities, namely, the axial anomaly and the definition of 75 
in n dimension |[Tl|. It is well known that the polarized triangle diagram for ^<K/q (see 
Fig. 1) has an axial anomaly. There are two extreme ultraviolet regulators of interest. One 
of them, which we refer to as the chiral-invariant (CI) factorization scheme, respects chiral 
symmetry and gauge invariance but not the axial anomaly. This corresponds to a direct 
brute-force cutoff ~ /i on the fcj_ integration in the triangle diagram ( i.e. fc^ < yU^) with fc_|_ 
being the quark transverse momentum perpendicular to the virtual photon direction. Since 
the gluonic anomaly is manifested at fc^^ — i> 00, it is evident that the spin-dependent quark 
distribution [i.e. Ag^^^(x)] in the CI factorization scheme is anomaly-free. Note that this is 



the fc^-factorization scheme employed in the usual improved parton model |13 



The other ultraviolet cutoff on the triangle diagram of Fig. 1, as employed in the ap- 
proach of the operator product expansion (OPE), is chosen to satisfy gauge symmetry and 
the gluonic anomaly. As a result, chiral symmetry is broken in this gauge-invariant (GI) 
factorization scheme and a sea polarization is perturbatively induced from hard gluons via 
the anomaly. This perturbative mechanism for sea quark polarization is independent of the 
light quark masses. A straightforward calculation gives |1^, |12[ 



A0;/a(^)Gi = ^<P\Ux) 



CI 
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X) 



(7) 



where the term — (1 — x) originates from the QCD anomaly arising from the region where 
k^^f 00. Two remarks are in order. First, this term was erroneously claimed in some early 



literature [^, |^ to be a soft term coming from fcj 



mz. 



(see below). Second, although the 



quark spin distribution inside the gluon A0^/J^(a;) cannot be reliably calculated by pertur- 
bative QCD, its difference in GI and CI schemes is trustworthy in QCD. Since the polarized 
valence quark distributions are /^^-factorization scheme independent, the total quark spin 
distributions in GI and CI schemes are related via Eq. (6) Hi 



Aq{x,Q^)Gi = Aq{x,Q^)ci 



as{Q' 



TX 



x)®AG{x,Q' 



For a derivation of this important result based on a different approach, namely, the nonlocal 
light-ray operator technique, see Miiller and Teryaev |jl^. The x dependence of the anomaly 
effect is thus fixed. 

The axial anomaly in the box diagram for polarized photon-gluon scattering also occurs 
at fc^ -^ oo, more precisely, at kj_ = [(1 — x)/Ax]Q'^ with a; ^ 0. It is natural to expect that 
the axial anomaly resides in the gluon coefficient function ACq in the CI scheme, whereas 
its effect in the GI scheme is shifted to the quark spin density. Since ACq (x) is the hard 
part of the polarized photon-gluon cross section, which is sometimes denoted by gi{x), the 
polarized structure function of the gluon target, we have 



AcS\x)=gf{x)-A<j>^;^^{x). 
It follows from Eqs. (^ and (9) that 

ACS\xhi = ACS\x)ci + -(1 - x). 



TT 



(9) 



(10) 



It has been argued that the GI scheme is pathologic and inappropriate [Q based on the 
observation that a direct calculation of gf{x) using the mass regulator, for example, for the 
infrared divergence gives 



^f(x)=^(2x-l)(ln% + ln^^ 



-1 +^(l-x) 



;ii) 



where the last term in Eq. (fiA] ) is an effect of chiral symmetry breaking and it arises from 



the soft region k'j_ ~ m?. By comparing (|TOD with ([TT|), it appears that ACq {x)gi, which 
is "hard" by definition, contains an unwanted "soft" term. This is actually not the case. 
Choosing a chiral-invariant cutoff on the fcx-integration, a perturbative QCD evaluation 



yields (see [Q for a review on the detail of derivation) 



.(1) 



AC'a''{x)ci 



27r 



{2x- l)(ln 



1 — X 



X 



(12) 



The ^(1 — x) term disappears in ACq {x)ci, as it should be, but it emerges again in the GI 
scheme due to the axial anomaly [see Eq. (|TOD or (|^)] and this time reappears in the hard 
region. That is, the gluon coefficient AC^j (a;)Gi is genuinely hard. 

It is easily seen that the ffist moments of ACg{x), J2q^(l{x) and gi{x) are given by 



and 



'dxAd^\xhi = 0, f'dxAd^\x)ci = -^, (13) 
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ASGi(g^) = ASci(g^) - ^^^^^AG(g^), (i4) 



p/ r^2\j -'- v^ 2 / A /'/n2\ '^sKSi 



Tl^j^gl{x^Q')dx = ^E^; Agci(g^)-^^AG(g^ 



= ^EejAgciig^), (15) 

where AG = /q AG{x)dx, and we have neglected contributions to f^f from A(j) L and AC^^\ 
Note that ASGi(g^) is equivalent to the singlet axial charge {p, s\J^\p, s) . The well-known 
results (12-14) indicate that Ti receives anomalous gluon contributions in the CI factorization 
scheme (e.g. the improved parton model), whereas hard gluons do not play any role in Fj* 
in the GI scheme such as the OPE approach. From (14) it is evident that the sea quark or 
anomalous gluon interpretation for the suppression of T^ observed experimentally is simply 
a matter of convention ||18[. From Eqs. (P) and (|10D we see that Aqf(x) + ACq (x) 0AG{x) 



and hence gi{x) is independent of the choice of the /c^-factorization scheme, as it should be. 
We would like to stress that physically, the GI and CI A;_|_-factorization schemes are exactly 
on the same footing, though philosophically one may argue that it is more natural to have 
ACg include all short- distance contributions. 

The MS scheme is the most common one chosen in the GI factorization convention. 
However, the quark coefficient function ACJj^\x) in the dimensional regularization scheme 
is subject to another ambiguity, namely, the definition of 75 in n dimension used to specify 
the ultraviolet cutoff on A0^A (see Fig. 1). For example, ACJ^^\x) calculated in the 75 
prescription of 't Hooft and Veltman, Breitenlohner and Maison (HVBM) is different from 



that computed in the dimension reduction scheme W^. The result 



AC«(x) = C«(x)-^(l + x) (16) 

usually seen in the literature is obtained in the HVBM scheme, where Cq{x) is the unpolarized 
quark coefficient function. Of course, the quantity Aq^^\x) + ACg{x) ® Aq^^\x) and hence 
gi{x) is independent of the definition of 75 in dimensional regularization. 
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In order to determine the Q^ evolution of the polarized structure function gi{x,Q'^) to 
NLO, it is necessary to know the two-loop splitting functions APj- (x) in the NLO evolution 
equation. Since the complete results for AP^^ [x) are available only in the MS scheme, a GI 
factorization scheme, it is natural to ask how do we analyze the DIS data of gi in the CI 
scheme ? One possibility is proposed in ||T6| that the evolution of the parton spin distributions 
Ag(a;, Q^)gi ^^iid AG'(x, Q"^) is first determined from the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution equation, and then Ag(a;,(5^)ci in the CI scheme is related to 
Ag(a;,Q^)Gi and AG'(x,(5^) via Eq. (||). The other equivalent possibility, as advocated 
by Miiller and Teryaev |]1^, is to make a renormalization group transformation of NLO 



splitting kernels and coefficient functions from the GI scheme to the CI one. In addition to 
the coefficient functions [cf. Eq. (10)] 

AC«(x)ci = AC«(a;)Gi, AC«(x)ci = AC«(a;)Gi - A{x), (17) 

where A{x) = —(1 — x), the NLO splitting functions in the CI scheme can be obtained by 
applying Eq. (|) to the spin- dependent evolution equations: 



|a,„.,M) ^ ^AP-' 



-AgNs(x,t) = ^AP„f ^x) ® AgNs(x,t), 



d_ f Ags(x, t)\aAtlf AP^s^(x) 2nfAP,G{x) \ / Aqs{x, t) \ 
dt[AG{x,t)) 2n [aPg,{x) APgg{x) J [AG{x,t))' 

where t = ln((5^/A^ ), and the indices S and NS denote singlet and non-singlet parton 
distributions, respectively. It is straightforward to show that 

AP,7(a;)ci = ApN/(x)Gi, APc;,(a;)ci = APG,(a;)Gi, 
AP^g{x)ci = API{x)gi + A{x)(E)APg,{x)gi, 

APgg{x)ci = APgg{x)gi - A{x) (^ APGg{x)Gi, 

a 3 
2nfAP,G{x)Gi = 2nfAP,Gix)Gi-^A{x) 

-A{x) ® [ApS - APgg + APg, ® A]{x)^^, 



(19) 



where 



"-("-5"^)^ £("-¥"')-■■■ P°) 

is the usual /^-function. The above results are first obtained by Miiller and Teryaev [|l^. In 
short, Eqs. (|l^) and (|19D provide the NLO coefficient and splitting functions necessary for 
the CI factorization scheme. It is easy to check that Agci = Jo dxAq{x)ci does not evolve 
with Q^ and that hard gluons contribute to Ff in an amount as shown in Eq. (14). 



Several Cl-like schemes were proposed in [^ in which the singlet anomalous dimension 
and the first moments of the gluon coefficient function are fixed to be 

Jo "- '"' 27r^ 



dxAcS\x) = -^, A7^;4^ = 0, (21) 



where 



A7I;. = /' AP,,(x)x"-Mx = A7i;^'" + ^A^ir + ■■■■ (22) 



ZTT 



The remaining moments of the coefficient functions and anomalous dimensions are then 
constructed by modifying the counterparts in the MS scheme. Three of such schemes, namely 
the Adler-Bardeen (AB) scheme, the off-shell scheme, and the Altarelli-Ross scheme were 
considered in 0. For example, higher moments in the AB scheme are fixed by requiring 
that the full scheme change from the MS scheme to the AB scheme be independent of x, so 
that the large and small x behavior of the coefficient functions is unchanged. Evidently, this 
is in contradiction to the (1 — x) behavior of the anomaly effect shown in Eq. ([l^); that is, 
though the first moments of Aq{x) and gi{x) in the AB scheme are in agreement with the CI 
factorization scheme, this is no longer true for other moments. In principle, a measurement 
of higher moments of spin-dependent parton distributions will enable us to discern between 
CI and AB schemes. Since the AB scheme is not consistent with perturbative QCD for 
the X dependence of the renormalization group transformation for coefficient and splitting 
functions, it is thus not pertinent to use this factorization convention to analyze the polarized 
structure functions in NLO. 

To conclude, since the x dependence of the axial-anomaly effect is fixed, the transfor- 
mation of spin-dependent coefficient and splitting functions from the MS scheme to the 
improved parton model (or chiral-invariant) factorization scheme in which the axial anomaly 
is shifted to the gluon coefficient function are uniquely determined by perturbative QCD 
for a given 75 prescription, contrary to the commonly used Adler-Bardeen scheme. We thus 
believe that the CI scheme should be used, instead of the less consistent and ambiguous AB 
scheme. In order to see the factorization scheme independence of Ag^,(x, Q^), AG{x,Q'^) 
and gi{x,Q'^) and demonstrate the effect of scheme dependence for Aqs{x,Q'^), AS(a;,Q^), 
it is urged that the QCD analysis of polarized structure functions in NLO should be carried 
out in both MS (or any gauge-invariant) and chiral-invariant factorization schemes. The Q^ 
evolution of parton spin distributions in the latter can be obtained by either studying the 
Q^ evolution first in the MS scheme and then applying Eq. (P) afterwards or solving the 
DGLAP equation directly in the CI scheme using the splitting functions given by (p!9|). 
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